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Interaction of 3 solitons for the GKdV-4 

equation 
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Abstract 

We describe an approach to construct multi-soliton asymptotic so¬ 
lutions for non-integrable equations. The general idea is realized in 
the case of three waves and for the KdV-type equation with nonlin¬ 
earity u‘^. A brief review of asymptotic methods as well as results of 
numerical simulation are included. 
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1 Introduction 

1.1 Statement of the problem 

It is well known that an arbitrary number of solitary waves collide for inte- 
grable non-linear equations in an enormous manner: they pass through each 
other almost as linear waves. The aim of this paper is the consideration; can 
we believe that such type of interaction is conserved (in a sense) for some 
essentially non-integrable equations? 

Needless to recall that the integrability implies both the possibility to 
find exact solutions of complicated structures and that the equation has 
some special properties. Conversely, the non-integrability implies that we 
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do not have, at least nowadays, neither explicit solutions, nor special useful 
properties of the problem. 

We consider the general techniques and averaging method by a specific 
but typical example of the Generalized Korteweg-de Vries-4 equation with 
small dispersion, that is: 


du du^ 2 * 9 ^'^ 
dt ^ dx ^ ^ dx'^ 


0, m = 4, X E t > 0, 


( 1 ) 


where e << 1 is a small parameter. 

The GKDV family ([T]) contains integrable (m = 2 and m = 3, that is the 
KdV and MKdV equations) and essentially non-integrable (m > 4) equa¬ 
tions. The last means that there is not any method at present to construct 
exact solutions of the Gauchy problem with more or less general initial data. 

More in detail, there is known that CD for m > 6 is unstable, whereas 
the case m = 5 is conditionally stable (unstable for solitons) [H I2]- As for 
m = 4, this case is stable and there are known some exact particular solutions 
including the solitary wave: 


u = Aijj (^-—, uj{r]) = ccosh (2) 


where /3 > 0 is an arbitrary number, c is such that 

poo 

/ u{ri)dr] = 1, 


( 3 ) 


and 

V = 

Here and in what follows we use the notation 





k > 1, a '2 



( 4 ) 

( 5 ) 


and the identities 



( 6 ) 


In view of the general wave propagation theory, there appear questions 
both about the stability of the solitary wave solution ([2]) with respect to 
small perturbations of the equation, and about the character of the solitary 
wave collision. 
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1.2 Prehistory: single-phase asymptotic solutions 

The non-integrability implies the use of asymptotic approaches. We consider 
waves of arbitrary amplitude (of the value 0(1)), treating the dispersion e 
as a small parameter. Thus, there appear ’’fast” x/e, t/e and ’’slow” x, t 
variables. Of course, it is possible to rescale 7] = x/e, r = t/e and pass to 
’’fast” T], T and ’’slow” erj, er variables. However we prefer the hrst version. 

The modern asymptotic technique, which is based on ideas by Poincare, 
van der Pol, Krylov-Bogoliubov and others, has been created hrstly by G.E. 
Kuzmak (ODE, 1959, [3]) and G.B. Whitham, (PDE, 1965, [HIS], see also 
[B]) for rapidly oscillating asymptotic solutions of non-linear equations. The 
famous Whitham method deals with a Lagrangian formulation and allows to 
hnd slowly varying amplitude and wave number of non-uniform wave trains. 
This approach determined the development of the non-linear perturbation 
theory in 1970s. At the same time, the passage from the original equation 
to the Lagrangian seemed to be artihcial. For this reason J.G. Luke (1966, 
0 ) created a version of the Whitham method, which allows to construct 
asymptotic solutions with arbitrary precision appealing directly to the origi¬ 
nal equation. More in detail, for the equation L(m, eut, eux ,...) = 0 we write 
the ansatz 

u = Yo{T,t,x)+eYi{T,t,x) + (7) 

where r = S{x,t)/e, Yk^r + T,t,x) = Yk{T,t,x), T = const, and S{x,t), 
Yk{T,t,x) are arbitrary functions from C°°. Since 

edtYk[S{x,t)/e,t,x) = {Std rYk{T,t,x) + edtYk{T,t,x)}\r=s/e, 

we obtain the chain of ordinary (with respect to r) equations, the hrst of 
them is non-linear, 

L{Yo{T,t,x),StdrYo{T,t,x),S^drYo{T,t,x),...) =0, (8) 

and others are non-homogenous linearization: 

L' {Yq{t, t, x), Std r, • • •) W(r, t, x) 

= Ffc(Yo('r,^,a:),... ,Tfc_i(r,t,x)), k>l. (9) 

It is assumed that (|H]) has a T-periodic solution. Then there appear the 
orthogonality conditions 

Fk{YQ{T,t,x),... ,Yk_i{T,t,x))ZidT = i = l,...,/, (10) 
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which guarantee the solvability of Q in the space of T-periodic smooth 
bounded functions. Here {Zi,i = 1,...,/} is the kernel of the operator 
adjoint to L'. Moreover, flTOj) allow to define the phase S{x,t) and all the 
’’constants” of integration of the equations ([8]), ([9]). 

It seemed that the same procedure can be used to construct a perturbed 
soliton-type solution (with trivial alterations). However, it is not true, and 
a mechanical repetition of the Whitham construction leads to some ’’para¬ 
doxes” and senseless solutions (see, for example, [8], pp. 303 - 306). The 
situation has been improved by V. Maslov and G. Omel’yanov (1981, [9], see 
also HD]). A little bit later a similar construction has been developed by I. 
Molotkov and S. Vakulenko (see e.g. HD). To illustrate the modification [9] 
let us consider the perturbed GKdV-4 equation ([T]), 


du du^ 2 
dt dx ^ dx^ 


( 11 ) 


where R = R{x, t, u, eu^, ■ ■ ■) is ’’small” in our scaling and R\u=o = 0. 

To find a self-similar soliton-type asymptotics we restrict the soliton part 
of the solution on the zero-level set of the phase S{x,t) = x — ip{t) + 0{{x — 
<p(t))^). This allows to avoid the appearance of some nonuniqueness effects 
(see (To], pp. 24 - 26). Next we take into account that a smooth small ’’tail” 
can appear after the soliton. Therefore, instead of ([^ we write the ansatz in 
the form: 

u = Yo{tR) + eYi{TR,x) + ..., ( 12 ) 


where r = {x — (p{t)) /£, Y^ are smooth bounded function such that Yq{t, t, x) 
tends to 0 as r ^ ±cxo, Ykir, t,x) —)■ 0 as r —)■ -|-cxo, and Tfc(r, t, x) —)■ YfT^x, t) 
as r —)■ —cxo for k > 1, and ip belongs to C°°. 

Similar to (IH|), substituting flT2ll into flTTll we obtain the nonlinear model 
equation 


dYo dY,^ d^Yp 
dr dr dr^ 

and define the shape of the leading term in ffT2|) . 


(13) 


Vo = A{t)u{(3T -F (pi{t)), /3^/^(t) = 'yA{t), 


(14) 


as well as the similar to (jlj) relation 

^ = “ 4 . 4 =>( 0 . 
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( 15 ) 







Here is a ’’constant” of integration. Next to find the deficient relation 
between ip and A we consider the hrst correction Yi freezed on the soliton 
front X = (fit). Denoting hi(r, t) = yi(r,t, x)\x=^p{t), we pass to the equation: 

d r d‘^Y 'i 

— I - (pA + AY,^Y, + Fo, Fo,, Yoxr, ...)-Yot. (16) 


Respectively, to guarantee the existence of the desired correction Fi, we ob¬ 
tain the following conditions: 


^Ottj 

.. .)dr. 

(17) 

Yqtt^ • • • ) 

- Yot^dr. 

(18) 


d 

dt 


Calculating the integrals in ffT7|) . we complete flT5|) by the equation 

d 

where 


(19) 


91 ^/ ij(r!)R(ip,t,A(j(r!),A/3ij(r!)^,Af3^aj(r!)^ri^...)dr!. ( 20 ) 


This allows us to determine the phase and amplitude dynamics. The equa¬ 
tions flTHl) . flT^ have been called ”Hugoniot-type conditions” [9] since they 
do not depend on e, whereas the solitary wave fflTl) Yq(^{x — ip{t))/e,t) dis¬ 
appears (in V sense) as e —)■ 0. Let us recall that the Rankine-Huginiot 
conditions remain the same both for parabolic regularization of shock waves, 
and for the limiting non-smooth solutions. 

Furthermore, returning to the asymptotic construction, we note that fllSI) 
implies the equality 

1 d A 

a^A^Y^ = - J R{ip,t,Auj{r]),Af3uj{r])r^,Af3^u{r])rjr„.. .)dr] - (21) 


where F^ = Fi|t-_,._oo- Now we integrate the equation flT^ and hnd the 
structure of the hrst freezed correction 


Fi(r,t) = y; (t)x(r,t) + Zi(r,t) -fci(t)FoV(r,t), (22) 
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where x and Zi are some fixed functions such that 
Zi —)■ 0 as r —±oo, 

X ^ 0 as r ^ +CX0, X ^ 1 as r ^ —cxo, 

and Cl is an arbitrary ’’constant” of integration. 

The next step of the construction is the extension of h"i(r, t) to hi(r, t, x) 
in the following manner: 


Fi(r,f,x) = (f,a;)x(r,t) + ^i(r,f) + Ci{t)YQ^{T,t), (23) 

where is a smooth function such that 

d u 

-^ = uiRu{x,t,0,..x<^p{t), t>0, (24) 

"^1 \x=ip{t) ~ ) t > 0. (25) 

Continuing the procedure we can easily construct the one-phase self-similar 
asymptotic solution with arbitrary precision. 

Let us note finally that self-similarity implies the special choice of the ini¬ 
tial data. In particular, the initial function Yi(r, 0, x) should be of the special 
form fl25D with arbitrary Ci(0) and arbitrary n^(a:, 0) under the condition 

ui{x,0)\x=ip{o) = ^'r(O). (26) 

If it is violated and, for example, u\t=o = A(0)a;(/3(a; — </?(0))/c), then the 
perturbed soliton generates a rapidly oscillating tail of the amplitude o(l) 
(the so called ’’radiation”) instead of the smooth tail eui{x,t) (see [12] for 
the perturbed KdV equation). However, eui{x,t) describes sufficiently well 
the tendency of the radiation amplitude behavior (see e.g. [I3|). 


1.3 History: two-phase asymptotic solution 

Concerning the solitary wave collision, this problem is much more compli¬ 
cated. Indeed, to describe the interaction of two waves of the form (I2|), we 
should look for the asymptotics as a two-phase function. 


u = W 


x-(pi 

e 


X-(P2 

e 



(27) 
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where hh(ri,r 2 ,t) has properties similar to the two-soliton solution of the 
KdV equation. However, to construct W we obtain a non-linear PDE, which 
is, in fact, equivalent to the original GKdV-4 equation ([1]). So, the existence 
of such asymptotics remains unknown. The same is true for any essentially 
non-integrable equation. Respectively, there is not any possibility to con¬ 
struct a classical asymptotic solution (that is, with the remainder in the 
G-sense). 

At the same time it is easy to note that the solitary wave solutions (soliton 
or kink type) tend to distributions as e —?■ 0. This allows to treat the equation 
in the weak sense and, respectively, look for singularities instead of regular 
functions. Obviously, non-integrability implies that we cannot hnd neither 
classical nor weak exact solutions. However, we can construct an asymptotic 
weak solution considering the smallness of the remainder in the weak sense. 

Originally, such idea had been suggested by V. Danilov and V. Shelkovich 
for shock wave type solutions (1997, [H]), and after that it has been devel¬ 
oped and adapted for many other problems (V. Danilov, G. Omel’yanov, V. 
Shelkovich, D. Mitrovic and others, [U] - [27] and references therein). We 
called this approach the ’’weak asymptotics method”. 

For the special case of soliton-type solutions we note now that they have 
the value 0(e) in the weak sense. Thus, the remainder for the leading term 
of the asymptotic solution should be O(e^) in the weak sense. However, the 
GKdV equations ([I]) degenerate to a first-order PDE in V' for this preci¬ 
sion. The same fact has been noted by Danilov, Omel’yanov, and Radkevich 
(1997, [28|) by the consideration of a free boundary problem. There has been 
suggested also a way about how to overcome this obstacle. Applying these 
ideas to the equation ([T|), we pass to the following definition of the weak 
asymptotic solution mi: 

Definition 1. A sequence u(t,x,e), belonging to C°°(0, T; C°°(M,[,)) for e = 
const > 0 and belonging to C(0, T; D'(M^)) uniformly in e > 0, is called a 
weak asymptotic mod Oxife"^) solution of (Qj) if the relations 

d r°° r°° All) 


d 

dt 
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u 


dfj , 

-T^dx -|- 3 
ox 




dfj 
d X 


dx = O(e^) 


(29) 


hold uniformly in t for any test function = 'iIj(x) G D(M^). 
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Here the right-hand sides are C°°-functions for e = const > 0 and piece- 
wise continuous functions uniformly in £ > 0. The estimates are understood 
in the C(0, T) sense: 

g{t,e) = 0{e^) t-)- niax \g{t,e)\ < ce^. 


Definition 2. A function v{t,x,e) is said to be of the value if the 

relation 

POO 

/ v(t,x,e)'ip{x)dx = 0{e^) 


holds uniformly in t for any test function G P(]R^). 


The sense of the relation fl28|) is obvious: it is the adaptation of the 
standard T>'-de£nition to asymptotic mod solution which belongs to 

C(0,T;Next we note again that (125]) cannot be a unique satisfac¬ 
tory condition since here has been lost the difference between the GKdV-4 
equation and the limiting hrst order equation (with e = 0). To involve the 
dispersion term into the consideration, we supplement fl28|) by the additional 
condition fl2^ . It can be treated as a version of fl28|) but for special test 
functions u'ijj{x), ^|J G which vary rapidly together with the solution. 

It is important also that fl22]) duplicates the orthogonality condition which 
appears for single-phase asymptotics. Indeed, the adaptation of the Defini¬ 
tion 1 to the perturbed equation (ITT]) implies the transformation of (l2^ to 
the following form: 


A 

dt 




,5-9 V'. 


U°^:f-dx -|- 3 
ox 


- 2 


' —OO 
POO 


du\ dip , 
Ox / Ox 


uRipdx = 0{e^). 


(30) 


Next for u of the form (1121) . (fTT|) we calculate the weak expansion: 



j^k 

u^{x,f)ip{x)dx = e— 



uj’'{r])ip{ip - 1 - e^)dx + 0{e^) 

Ak 

= eak—ip{^) + 0{e^), 


where k > 1. Thus 


u 


k 


j^k 

eak—S{x - 




(31) 


In the same manner we obtain 


(32) 

(33) 


= ea'2(3A^5{x — (^) + 

uR{x, t, u, eux, e^Uxx, • • •) = £^^5{x - ^) + Ov'{e^), 

where has been defined in fl20l) . Snbstitntion of flMl) - (l33|) into (l30|) implies 
the relation 

r dt^ 8 ^ ! n A2^ a/ \ 

n ■ * + “=5y ■ p(x-v>) 

Since 5{x — ^p) and 5'{x — (p) are linearly independent, their coefficients in 
fIMll should be equal to zero. Taking into account the identities (jlj) we obtain 
again the equations ffT5l) . ffT^ for the one-phase asymptotics ffT^ . 

Let us revert to the two-wave interaction. Following ca (see also [T8]L 
we present the ansatz as the sum of two distorted solitons ([2]), that is: 

M = ^ GiUJ > (35) 


where 

Gi = Ai + Si{T), ipi = ipio{t) + eipii{T), T = l3i{ip2o{t) - ipiQ{t))/e, {m) 

Ai are the original amplitudes and </?jo = + describe the trajectories of 

the non-interacting waves ([2]), (3i = {'jAi^A^ We assume that Ai < A 2 and 
Xio >> X 20 , therefore, the trajectories x = </?io and x = (p 2 o intersect at a 
point (x*,t*). Next we define the ’’fast time” r to characterize the distance 
between the trajectories (pio and we assume that 5'i(r), p>ii(T) are such that 

Si —>■ 0 as T —> ±cx), (37) 

(Pii —>■ 0 as T —>■ — 00 , tpii = constj as r ^ -l-oo. (38) 

It is obvious that the existence of the weak asymptotics (15^ with the prop¬ 
erties fl37)) . fl38|) implies that the solitary waves ([2]) interact like the KdV 
solitons at least in the leading term. 
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To construct the asymptotics we should calculate again the weak expan¬ 
sions for and {eu^Y. It is easy to check that 


At the same time 


^ G 

u = e'^ - ^i) + Ov'YY- 

i=i 


G 


(39) 


u‘^{x,t)'ilj{x)dx = I )dx (40) 

> 2=1 J—OO Pi 


2 G 1 G 2 I UJ ( /3i ^ ) CG ( (32^—-^ ) Y{x)dx. 


(41) 


We take into account that the integrand in fITT]) vanishes exponentially fast 
as |(pi — (^ 2 ! grows, thus, the main contribution gives the point x*. We write 


ipio = X* + Vi{t — t*) = X* + e 


Vi 


A(^2-W 


-r and ipi = x* + exu (42) 


where Xi = ViTj{p\{y 2 — hi)) -f ipi\. Next we transform the integral in (l4Ti) 
to the following form: 


£ I ^ Tj 

— / a;( 6 >i 2?7 - ai 2 )u{r])'ijj{x* 4- £^2 + )dr], 

P2 J -00 P 2 


(43) 


where 6*12 = (3i/(32, o'i 2 = f3i{^i — ^ 2 )/^- If remains to apply the formula 

f{T)S{x - ipi) = f{T)S{x - X*) - exif{T)6'{x - x*) + (44) 

which holds for each of the form fl42p with slowly increasing Xi 
/(r) from the Schwartz space. Moreover, the second term in fl44)) is Opi{e). 
Thus, under the assumptions fl371) we can modify fl39|) - flTT]) to the final form: 

“ ^ ^ ^ ~ ~ ~ ^Xi3'{x - X*)} -4 Ov>{eY, (45) 

i=i i=i 


2 2 

= ea2 ^ ~ Pi) + £ 02 ! ^ ^{2AiSi + S} 

i=l i=l 
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(46) 


G,G 


l'-^2 


(92 


A 2 ,i(cri 2 ) - X*) + Ov'ie ), 


where the convolution A 2 ,i(o'i 2 ) describes the product of two waves. In view 
of further applications we present such type of convolutions in the general 
version: 




r]^u'^ ’^{r]in)u}’'{r])dr], r]in = Oinf] - ain, Oin = ^, (47) 

Pn 


where 1 < < m, m > 2, j = 0 or j = 1, and we write Xm,k{o'in) 

simplifying the notation. 

To calculate the time-derivative of u with the accuracy it is 

enough to use the expansion fl45D . the assumptions fl371) . (I38|) . and to ap¬ 
ply the formula fl44D again. Thus, 


du 

Jt 


i=i i=i 

,2 , 

+ xXa‘js'(x - x") + Oxx(£^), 


(48) 


where i/^o = A (^2 - A)- 

On the contrary, to find d {v?‘)/dt with the same accuracy we should add 
to the leading term (1461) the next correction; 


Xi, 


- e^a2\ ^ ^{2AiSi + Sf 


i=l 


+ 2 —^-^X 2 A 2 ,1(0-12) + ^A2j(cri2) j jd'{x — X*) + Ox'/(£^). 

Now we obtain 

= a2^|2^GiG2A2,i(cri2) ^ ^( 2 AA Sf)^ 6 {x - x*) 

~ ^ 2^GiG2(x2A2,1(0-12) + ^A2^i(0-12)) 


(49) 


(50) 
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'^Xii^ASi + S"^) \5\x - X*) - ea2 ^ Vi^6'{x - ipi) + Ov'{e^)- 


A? 


2=1 


A 


2=1 


A 


Calculating weak expansions for the other terms from the left-hand sides of 
(EH]), ll2^ and substituting them into fl28l) . fl2^ we obtain linear combinations 
of 5'{x — ipi), i = 1, 2, 5{x — X*), and 5\x — x*). Therefore, we pass to the 
following system of 8 equations for 8 unknowns: 


F,,AA,A,A) = 0, j = 1,2, * = 1,2, (51) 

^g,(A,^2,cri2) = 0, j = l,2, (52) 

^ = F,(A,^2,cri2), J = l,2. (53) 

The hrst four algebraic equations fl5Tll imply again the relations (jlj) between 
Ai, (]i, and Vi. Furthermore, functional equations fl52|) allow to dehne Si with 
the property fl37|) . whereas an analysis of the ODE fl53|) justihes the existence 
of the required phase corrections ipn with the property fl38|) . By analogy with 
(1T5|) . (IT^ we call (|^ - (j53l) the Hugoniot-type conditions again. Results of 



.r 


Fig. 1: Evolution of two solitary waves for e = 0.1 

numerical simulations [29l [30] conhrm the traced asymptotic analysis (see 
Figure 1). 

Finally let us note that the two-wave problem for the equation ([T]) has 
been considered recently in the framework of another approach but for small 
amplitudes [3T] . 
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1.4 The next problem: three wave interaction 

Further numerical investigation of the GKdV-4 equation showed that N soli¬ 
tary waves collide elastically (in the leading term) at least for N < 5, see 
[29l [30] and Figures 2-4. Thus, there appears the problem of describing the 



Fig. 2: Evolution of the soliton triplet for e = 0.1 


interaction of three (and more) solitary waves. It seemed that it was enough 
to repeat the same procedure as above, now for three waves. However, it is 
easy to recognize that the corresponding 4 equations (|52|) . (]53|) will contain 
now 6 free functions. Obviously, this can not be any adequate description of 
the solution. 

Therefore, we should transform the conception of the asymptotic solution. 
To do it let us recall two hrst conservation laws for the equation ([T]) (in the 
differential form): 



(64) 



( 55 ) 


Comparing the left-hand sides of fl5T|) . fl5^ with fl28|) . fl29|) we conclude that 
Dehnition 1 calls a function to be a ’’weak asymptotic solution” if it satishes 
the conservation laws (j54j) . (l55l) in the sense Od'{s^)- Next note that pertur¬ 
bations of ([I]) imply corresponding perturbations of the conservation laws. 
For (ITTl) . instead of fl5T|) . (I55D . we have the relations 



( 56 ) 
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X 


Fig. 3: Evolution of 4 solitons for e = 0.1 


du'^ 


A 

dx 


-u 


5 




+ e" 




■1 


= 2Ru. 


(57) 


\dx/ ' dx"^ 

Reverting to the single-phase asymptotic solution (|T2|) one can easily es¬ 
tablish that the orthogonality condition flTT)) . and therefore the equation (IT^ . 
is the integral form of flFT)) . calculated for flT^ with the accuracy Ov'{£^)- 
At the same time flTSl) . and thus the equality fl2T]) . is the integral form of 


mod calculated for flT^ . where Yi has the form fl23ll and Ui 

satishes the equation (EH). Note also that the second ’’conservation law” 



Fig. 4: Evolution of 5 solitons for e = 0.1 

fl57)) has been used in the one-phase situation to dehne the principal term of 
the asymptotic solution, whereas (l56ll has been used only to dehne the hrst 
correction Yi. 
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Therefore, we see that to dehne the principal asymptotics term there has 
been used only one conservation law for the single-phase solution, and two 
conservation laws for the two-phase solution. So it is natural to assume that 
to construct a three-phase asymptotics we should add to fl5T|) . fl55|) the third 
conservation law, namely 


¥i 


/ 

\ dx) 



d 


dx 


,dudu 


— + £ 


dt dx 


, d'^u \' 
dx'^J 


0. (58) 


2 Asymptotic construction 


Let us consider the equation ([T|) with the Cauchy data 

3 


u\t=o = ^AiUi ( A 
i=l 


X - Xio 

e 


(59) 


where Ai < A 2 < A^, xio » X 20 » x^q. 

The arguments considered in the previous subsection imply the following 

Definition 3. A sequence u(t,x,e), belonging to C°°(0,T;C°°(M,),)) for e = 
const > 0 and belonging to C(0, T; T>'(M,^)) uniformly in e, is called the weak 
asymptotic mod Ox>'(s^) solution of the problem lljf), / f5PI) if the relations 
and 


¥i 

+ 


2e^ 




,dudu 
dt dx 


u 


(60) 


+ + U 


.d‘^z 


dx'^ 




dx'^) dx 


dx = 0(£^) 


hold uniformly in t for any test function f/ = 'ijj{x) G P(M^). 

To construct the asymptotic solution we write the ansatz in the form 
similar to fl35ll . namely 


u = 




i=l 


X - ipi 


(61) 


where the same notation and hypothesis fl36ll - fl38l) are assumed with obvious 
corrections: the ’’fast time” is dehned now using the distance between the 
hrst and third trajectories. 


= A(</?3o(i) - 


(62) 


15 







and we suppose the intersection of all trajectories x = (fioit), i = 1,2,3, at 
the same point {x*,t*). 

The technic of our approach has been explained in Subsection 1.3. So we 
clarify here some new detail only using as an example All others explicit 
formulas for the asymptotic expansions are presented in the Appendix. 

Lemma 1. Let the assumptions [3l\ ), foru of the form be satisfied. 
Then the following asymptotic expansions hold: 

3 

- ipf) + ^md{x - a:*)| + (63) 

i=l 

= amWQ—^o[x -X ) - earn 2_^ ViKlfy 'd {x - ipi) 

i=l 

d ^ 

- — ^ -X*) + Ovfe'^), (64) 


where m > 1, 


Am 

i’o = A(n - ro, y”’ = K\f = K\p = /fi”' - yy. 

3 

A'jl ^ Rm,ln + -Rm,123) 

2=1 l,n 

3 

^ + '^^{XnRm,ln + Cm,In) + XsRm,123 + C'm,123) 


2=1 l,n 

m—1 


Rm,ln = A 


k=l 


m—1 j — 1 


Rm, 123 = AA J2J2^'rnC^^r K2^3>^^ 
j=2 k=l 

m—1 

c„j„ = Y1 


k=l 


m—1j—1 


Cm,123 — AA 3 

j=2 k=l 


(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 

(71) 
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are the binomial eoefficients, the notation (f^TP has been used, and 


fin = fl2 + fl3 + /23, cr^n = f^l— - Xi = + ifn. (72) 

Furthermore, 

1 

^m,i 23 ^ = — / v"^'^~^{Vi3W~’"{V23)uJ^{v)dv, * = 0 or i = l. (73) 

J—oo 

To prove the lemma let us separate all the terms of u™ into three groups: 
one-phase, two-phase and three-phase functions: 

3 m—1 m—l j—1 

+ ctyr'‘yt + T.T.cicfYr’rtx. ( 74 ) 

2=1 Z,n k=l j=l k=l 

where Yi = Giu{^j3i {x — ipi)/e). Now considering n™ in the weak sense we 
change the variable: x = Lpi + erj/(di, x = ipn + ^'h/t^n, and x = ips + erj/fd^ re¬ 
spectively for the integrals of the groups. Next, preparing the same transfor¬ 
mations as in Subsection 1.3 and applying fHT)) again we pass to the formula 

dsnu- 

In the same manner one can prove the following proposition: 

Lemma 2. Let the assumptions (PT7| ), (E^) be satisfied for u of the form 
Then the following asymptotic expansions hold: 

3 

{euxf = - ipi) + 91(i),2(5(a: - a;*)| -h (75) 

2=1 

^(£ m ^)2 = - X*) - ea2^fiiViKl^'’d'{x - ipi) 

2=1 

- ea'2fio-^[^l3iF?Tii + + Ovfe^), (76) 

2 = 1 

3 

(e^Mxx)^ = - ipi) + D\^2),2d{x - a:*)| Od'(£^), (77) 

2=1 

3 

= £ 023 ! ~ - (fii) + 2.6{x - x*)| -h Ovfe'^), (78) 
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e'^u^ut = -ea^^ 6{x - (pi) - - x*) + Ot>'(£^), (79) 

2=1 

where a' 2 , ipo, defined in (E^, ^ = '>po{& + &g) + 937, &, 

&G, 93t, and other notation are deciphered in Attachment, Subsections 6.1 
and 6.2. 

Now we substitute the expansions ([63]), fl64|) . and fl75l) - d79D into fl28|) . (1^ . 
(1601) and obtain the similar fl5T]l - fl5^ system. Namely, the algebraic system 


for each i = 1, 2, 3: 

- KA'S’ + atK'i = 0, (80) 

- a.y.Kf + - SaiftUl'® = 0, (81) 

- a'yyy'i + josK aS> + 24ftV./s:® 

- ogA-® + 8o23/?f Ag> - <,;/?*Ag> = 0, (82) 

the system of functional equations: 

Z ^7 = 0, (83) 

2=1 

972 = 0, (84) 

0297(1),2 - ^05975 = 0, (85) 

and the system of ordinary differential equations: 
d ^ 

- H + “4^4 = 0, (86) 

2=1 

7 3 

- + 97^'^} + ^05975 - 3a'97(i),2 = 0, (87) 

2 = 1 

- «2 ( Vil + 97|;5 2 ) + Tsrfo Vn + 97^^^) } 

2=1 2=1 

+ — 0 - 8^8 ~ 2a23-C — 02 97(2),2 = 0. (88) 


Let US overcome the hrst obstacle: for each i the system (l80j) - (j8^ of three 
equation contains only two free parameters Ai, Vi. 
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Lemma 3. Let uj{r]), Ai = A{(3i), and Vi = V{(di) be of the form 
Then the equalities are satisfied uniformly in (3i > 0. 


Proof. Obviously, equations (ISO]) . (I5TD coincide with flHT]) and imply again 
the formulas ([2])-([2]). Substituting them into we transform it to the 

following form: 

0204 - as - a'a'y® + 80237 ^ = 0. (89) 

Next we note that a;( 77 ) satishes the model equation 


dr]'^ 


= 040; — 


(90) 


Multiplying fl90|) for u” and integrating, we obtain the identity 


4023 = 7^02 + 0402- 


(91) 


On the other hand, integrating the squares of the left-hand and right-hand 
parts of fl90|) . we pass to another identity: 


as = 7^02 — 0204 -|- 20405 . (92) 

This and ([ 6 ]) verify the equality flH^ . □ 

Since the system of six equations fl83l) - fl88ll contains six free functions, we 
obtain the hrst formal result 


Theorem 1. Let the system have a solution which satisfies the 

assumptions of the form |T7| j, (E^j- Then the solitary waves m collide 
preserving mod Ox>'{£^) the KdV-type scenario of interaction. 


Moreover, similar to the Rankine-Hugoniot condition, which is simply the 
conservation law for the shock-wave solution, the Hugoniot-type conditions 
fl80|) - fl88|) imply the verihcation of some conservation laws: 

Theorem 2. Let the assumptions of Theorem 1 he satisfied. Then the ansatz 
/ fhi]) is a mod asymptotic solution of the equation (Q]) if and only 

if (E2P satisfies the conservation laws 


fL 

dt 





2 

^ dx) 



dx 


0, (93) 
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and the energy relations 
d 


dt 

d 


xudx — / u dx = 0, 

J —OO 


' —OO 

/ OO o poo 

xu^dx -' 

5 


u^dx + ?> I (] = 0, 


— OO J —OO 

^{rx(e^^)^dx-^-rxu^dx 

dt \J_^ \ dx) 5 J_^ 

^ 2 r dudu^ r ( 4 od^u^^ 
+ 2 ^ / iu^ + e 

J_^ dt dx J_^ 


dx"^ 


^ ' dx = t). 


To prove this conclusion it is enough to rewrite the equations 
in the integral form. 


(94) 


3 Analysis of the Hugoniot-type conditions 

3.1 Transformat ions 

Let us pay the attention to the equations (l83l) - (|85l) . Normalization 

implies 


n^tn/ o—i. 

= where A. = 1 + 9‘f 


We denote ^ and obtain: 


(95) 


(96) 


= Y. - 1 ) + I] 


2=1 

m —1 


L^n 


E fyk n 
k=l 

m —1 j —1 


k ^ 2 (m—fc)/ 3 ^ 2 fc/ 3—1 * m—k \ k 


9S''‘^‘Ar‘’A;A„^,(<T,„), 


j =2 k=l 


kf,‘ 2 im-j)/ 3 Q 2 {j-k )/3 jym-j jyj-k jyk ^( 0 ),{j,k) 
13 23 1 2 3 ?t?/j 123 


(97) 

(98) 

(99) 
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This and similar formulas for ^{i ),2 (see Attachment) allow 

us to transform fl83|) - fl85|) to the following form: 


3 


^ Ki = 0, 

2=1 


(100) 

ef(X - 1) + 2 A, A„ X2.M„) 

2=1 l,n 

= 0, 

(101) 

E- 1) + 2EAi A„A'",Vin) 

2=1 Ln 

4— 

- 3^H5 = 0, 

(102) 


where the equalities ([6]), the notation (I97j)-(l99]), and (I140p have been used. 

Next let us simplify the equations (|86l) - (l88|) . We note hrstly that in view 
of (183|) and the identity 

f 3 l{xi - Xn) = (Tin (103) 

it is possible to eliminate Xi from the left-hand side of fl86l) . since 


i=l 


A'io'hii + XiAf’ = ^ + 


2=1 


^12 ^( 1 ) 


f32 '' ■ 


In the same manner, applying fl84)) and fl85|) . we simplify the equations (IHT)) . 
). Thus, we transform fl86p - (l55]) to the following form: 




2=1 

3 


2 = 1 




d 


L,n 

4 


23 


= F, 


E )v>ii + is J - 26 = 5, 


2=1 


(104) 


(105) 

(106) 


where 




f = a,D\,, F = ^%),2, ^ = 2m-^D^s-2^Q-^^Ri2),2, (107) 

02 On tto On 
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and the function ^ is described in Attachment (see formula 01421) b 

The second step is the elimination of (pn from the model system. To do 
it we divide ain into the growing (Win) and the bounded (Win) parts: 


and rewrite the identity 01O3j) : 


Thus 


<^12 


def /i/ / T /' T / 1 

= -^(Vi - Vn)T 

Wo 

(108) 

Pnl ) • 

(109) 

W 23 , 

.1 —-^ + P21- 

P2 

(110) 


Substituting OllOp into 01O4p we obtain 

+ pixP - + 2 


dr 


n dr /3i '' 


P 2 /?2 

Here and in what follows we use the notation 

3 

kIq for j = 1 and j = 2 . 


= 


+ ^. (Ill) 

n 


( 112 ) 


2=1 


Next we use the equalities 01101) . 01111) . and 

dl3 = 5 ‘i2 + 6*i25'23, (713 = (Ti2 + ^12<723) 


(113) 


and rewrite OlOSp . OlOOp as equations for new unknowns 5 'i 2 , 5‘23. After 
normalization fl9^ we pass to the following model equations: 


■00 


d 

dr 


(7i2 , (7i2 (723 (^23 

Pio -p, —1“ Pii “5-Pso -p -P 31 -p- 

Pl Pi P2 P2 


yic2,,„}=f’--/. (114) 

“ J n 

L,n 


■00 < 


A 

dr 


(7i2 , (7i2 


+ 2T 


d 


dr I 


lib!!) 
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(7i2 


A 

dr 


(723 , (723 




“ 2ri6G 


(115) 


where = ri^ + ( 2 T + (lio’)f the coefficients pf), cm, T are 
presented in Attachment (see formulas 0144p - 0148p ). 
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3.2 Asymptotic analysis 

To simplify the further analysis let us assume that 

^ 23 ^ = hi ^ 12 ^ = where a G [0,1) and /i is sufficiently small. (116) 
We look for the asymptotic solution of the system fll 00 l) - fll 02 p in the form: 

1 1 

Ki = k,2 = (117) 

where Xq and yi are free functions. Then fllOOp is satished, whereas fllOip 
and fllU 2 l) imply the system: 

2xo - Ai“|/i|l + /iA 2 , 1 ( 0 - 23 ) - /^"(l - A2,i(cri2) + ^yi) - /i^^^+“^/^A2,i(cri2)| 
= —2A2,1(0-23) — 2/i"A2,1(0-12) — 2/i^“''"A2,1(0-13) + 0 \{ y ^), ( 118 ) 

I? + 40/x^A5^4(o-23) + —y^xo^xo — 5y^~^°‘yi\5^4{a23) 

= —10A5^4(o- 23) — 10/i^’''“A5_4(o-i3) + (119) 

Here and in what follows we denote 

f{a,y) = Ox{y^) if max |/(o-,/i)| < c/x^ ( 120 ) 

<7 

and f{a, •) belongs to the Schwartz space. 

It is easy to see that the compatibility of the equations (IllSp and (11191) 
requires the condition: IOA 5 , 4 ( 0 - 23 ) = 7 A 2 , 1 ( 0 - 23 ) + Oa(/x“). 

Lemma 4. Let oj{ri) be of the form Then 

10A5,4(ain) = 7A2,i(az„) + 39in\fh^in) (121) 

for all indices /, n. 

To prove the lemma it is enough to use again the equation fl 6 SP and the 
identities ([5]), ([ 6 ]), fl67l) . 

Now we set 

Xq =—A2,1(0-23) +/x“a;i (122) 

and transform flllSp . flll9p to the hnal form: 

2x1 — ri 2 X/i =— 2 A 2 , 1 ( 0 - 12 ) — 2 /XA 2 , 1 ( 0 - 13 ) + Oa(/x^ “), (123) 
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^21X1 — ^/^A5,4(c’‘ 23)|/1 — —/^A2,l(cri3) — -/i^ °‘^n{o'23) + Ox{f/), (124) 

where 

ri2 = 1 + /iA2,l(cr23) - (l - A 2 , 1 ( 0 - 12 ) + ^yi) + 

IQ 

?^2i = 1 + ^/i^Aa,1(0-23) + OA(/i^'^“). 

Solving this system we obtain the asymptotic representation: 

xi = —/^(Aa,1(0-13) — A 2 ,i(o-i 2 )A 2 ,1(0-23)) + ( 125 ) 

yi = 2A2,1(0-12) (^1 + yU^+“(l - Aa,1(0-12)) j + OA(/i^“")- ( 126 ) 

Combining (11171) . (11221) . (11251) . and (I126p we conclnde: 

Lemma 5. Let there exist functions ipn, i = 1,2,3, with the properties / f5^) 
and let the condition Mldt) he realized. Then the system U00\) - M02i) has 
the unique solution 

1^1 = h'^A 2 ,i(o-i 2 ){l + — Aa, 1 ( 0 - 12 )) } + (127) 

^2 = ~h”A 2 ,i(o-i 2 ){l + /i^'''"(l — A 2 , 1 ( 0 - 12 )) } + 0 \{p?), (128) 

^3 = “h^Aa, 1 ( 0 - 23 ) + (129) 

such that Si = satisfy the assumptions [3l\ ). 

To complete the analysis we shonld prove the solvability of the system 
(11141) . (I115D . Taking into acconnt (195|) . (I116p . and (11271) - (11281) . we obtain: 

^ 11 ^-^ 12 ^ 12 ^=^!, (130) 

dr dr 

~ d(J^o ^ ~ d(7o‘i ~ 

E21—\ -^12.^22—1- = T 2 , (131) 

dr dr 

where the coefficients Eij and right-hand sides Fi are demonstrated in At¬ 
tachment, Snbsection 6.4. 

It is easy to calculate that 

det{Eij) = 6 ' 12 /nA, (132) 
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where 


A = ^ + + 0;,(/i(3+“)/2 ^ 

o 

Thus, we transform the system fll30p . fll3ip to the standard form 

= Mi 2 (r, (712, CT 23 ,/i)/A, = M23(r, (712, CT 23 ,/i)/A, (133) 

where 

Mi2 = ^(ct) = o'-^2,i(o'); (134) 

11^23 = ~2/i “A2 ,i(o'23) — -^^(< 712) + OA(/i^^ + /i“), (135) 

and the equalities ([S]), (1^ . (1^ . as well as the functional relation fll2ip and 
a&\7{(Tin) = a2a\\2,i{ain) - 7^023A4,3(2)(uzn) + {(Tin) (136) 

have been taken into account. 

According to the notation fllOSp and the first assumption of the form fl38l) 
we add to fll33p the ’’initial” condition: 

5‘i2|^^_^^0, 5-23|^^_^^0. (137) 

Since Mij vanish with an exponential rate as r — ?■ ± cxd , it is easy to prove 
the solvability of the problem fll33p . fll37p . Next we note that A 2 ,i(crj„) = 
A 2 ,i(^Zn' 7 ‘ + o'ln)- Since ^23 = we find from (11331) . (I135p that 

o' 23 {t) = for sufficiently large r, however it tends to the limiting 

value sufficiently slowly, with an exponent 0(/i^^'''“^/^). Conversely, taking 
into account that ai 2 = 0(/i®), we obtain that di 2 (r) = 0(1) for sufficiently 
large r and tends to the limit with an exponent 0(1). 

The last step of the construction is the return to the phase corrections 
(fii- In view of (IllOjl . (11111) it is obvious that the last assumption of the form 
(l38P is justified. This implies our main proposition 

Theorem 3. Under the assumption Ulf^) the asymptotic solution (El]) de¬ 
scribes mod Ot>'{£^) the KdV-type scenario of the solitary waves interac¬ 
tion. 
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4 Conclusion 


We looked for an approach to describe solitary wave collisions avoiding the 
use of explicit multi-soliton formulas. Surprisingly, we came back to the an¬ 
cient Whitham’s idea to construct asymptotics with the help of conservation 
laws and a reasonable ansatz, but in the framework of the weak asymp¬ 
totics method. In our case three conservation laws for three waves have been 
utilized. It is clear now how to generalize the approach: for N waves N con¬ 
servation laws should be used. On contrary, the existence of N conservation 
laws does not imply the existence of iV-soliton type solution since some very 
astonishing additional conditions appear to guarantee both the solvability of 
model equations (like 01211) 1 and the regularity of the solutions (like 0136p h 
Furthermore, some questions remain open, the hrst of them: how to choose 
the collection of conservation laws to describe Wsoliton interaction and is 
it possible to change conservation laws to reasonable energy relations? At 
the same time we can formulate the main result of the paper: there is not a 
sharp frontier between integrable and nonintegrable equations: similar sce¬ 
narios of the soliton interaction are realized, but with small corrections in 
the nonintegrable case. 
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6 Attachment 

6.1 Formulas for {£Ux)‘^ and (e^u^x)^ 


3 



3 



T]^u'{T]in)u'{T])dT], (140) 


(139) 



(141) 


and U)'{ri) = du{ri)/dri, uj''{r}) = d^ujirf)/dtf. 
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6.2 Formulas for e^u^ut and 




(2) dif^ 


i=l 

3 


dr 

'( 2 ) 


difii ^ difnl 


l^n 

^ 9?t;n = + Fn)Ai‘;;(azO, 


2 = 1 


Ln 


&G = 


u 

l^n 
4 


^ dGn „ dGi 

—;-(j'n 


dr 


dr 


. •a=-4 5;/3?Al? + 5;i;,„ + o, 


2=1 


l,n 


fi.„ = 5^ Cifl/?;/i:,'*-■’■>/(•«+'> A4j,2)(a,„), 

j=l j =0 

a = Y, Cl{l3il}2l3lK't’'’K2'K!i''Km 


i=i 


+ l3il3lKr^'K'^'K!,^>\n2) + l3!l32KrK2K3*Kj(3) 

GiG^PiP2Kt^^ Krh<rt 


j =2 k=l 


m=l 


Ao 7(cTzn) = — / u{r]in)u'{r])dr], 033= / a;'^(r7) (a;'(?7)) drj, 

^2 J —CO J —00 


'^4,j(l) (^/n) 

A4J(2)(0'in) = 
A5J(1) 

A5j(2)) = 
A5j(3) = 


023 J- 

O 23 J- 

2 


00 

00 


O 23 J- 

1 


00 

00 


O 23 J-00 

2 


A, 


023 J- 

1 


{Vin)uJ^ {ri)uj'\r]in)dr], 

uj^~^{Vin)(^\vW\v)dv, 
(^'^~^iVi3)^^iV23W'{v)dv, 
uj‘^~^{Vi3WiV23Wiv)dv, 
(^''iVi3)(^'^~^iV23)(^^iv)dv, 


00 

00 


6Jkm 


^{^13)^^ ^{V23)^^''{vW{Vm3)dr]- 


O23 J-c 
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6.3 Normalization 


«(i ),2 = E - 1) + 2 E A„A<“|((r,„), 

i=l l^n 

^ = fin - f (41’ + 03) + <,<r> = /3?A-2> + (142) 

= E (C 2 I!. - fc,M) - (143) 

l^n 

2.4 

<33 = ^ “ gT^-Rsjsj - 27^-^5,123, (144) 

i=i 

p,, = 3Rf^ = 0 if* ^ 3, M 1, (145) 

i=i 

eio = -nq^ + k\^ Y. - 2^1 ^ + ^i3 + 0) + 2/4'^ 4/, (146) 

i=i 

eii= r,(«‘;>+Q.)4-Al;’E«fo’' ^ = E 7A.?’4-2 E 0”. (l^^) 

j=l 1=1 In 

Cl=Cl2, C3 = C23, Ol = 0. (148) 

6.4 Asymptotic analysis 

^11 = -f2 + - r2/i='('+“)/2z'(ai2) - 2 ri/i''('+“)/ 2 ^'1(^12) 

+ 2fi/X^’^"^((TA2l(a + 012O'23))^1^^^^2 +Oa(4’^^“), 

Ei 2 = ri — r20il^^ — 2rifi^~^°'(^X2i{cri2)z\a23) — Ai2;'(cri3)) 

+ r2/i(^+^“)/^z'(a23) + 

i=l i=l 

^ - 4/i^A2l(cr23) + Oa(/^^'^^"), 

E 22 = -^ri + I^(3+i^)/2 + 4/i2(A2i(a23) - ^^'(^ 23 )) + C>a(/x^+“), 

4 = -^4{A2i(a23)(l +4+“)/^) +4+“)/^(2A43(a2^ 
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— A2i(cT23)) + OA(/i^)|, F2 — 2/i^|A2l(cT23) - (A 2 I (o'ls) 

A2i(c’‘12)A2i(0'23)) + (A43(c’‘23) ~ A 2 l(cr 23 )) + Oa(/^^)|- 
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